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The mobility edge is extracted from a non-perturbative analysis of F. Wegner's real matrix en- 
semble (RME) , iV-orbital model of electrons with broken time-reversal invariance moving in random 
potential. The replicon fluctuations around the zero-dimensional (0D), single metallic grain sad- 
dlepoint, are shown to destabilize diffusive regime beyond the low energy limit yielding a precise 
criteria of localization in terms of a single scaling parameter. 



Introduction- The mobility edge problem for disor- 
dered electron systems has an extensive history dat- 
ing back to the seminal work of Anderson in 1958, 
who designated a criteria of localization for the insu- 
lating phase based on the existence of discrete levels in 
the energy spectrum that are accompanied by normaliz- 
ablc cigenfunctions 1]. Localization signatures are found 
within the squared, impurity averaged Green's function, 
Una \G(r, r'; i)| 2 1], or more generally, higher moments 

of the wave function which also contain distinct non- 
analyticitiesQ- Although there are a variety of perturba- 
tive methods for disordered systems, all seem to break- 
down near the transition point (aside from the Bethe 
lattice), and the exact mobility edge trajectory has not 
yet been determined. 

An effective field theoretic approach was introduced 
by F. Wegner and collaborators !3| , that built on ear- 
lier ideas [J] developed from a replica field theory. Weg- 
ner identified a continuous symmetry of the iV-orbital 
version of Anderson's tight-binding model and a for- 
mulated the mobility edge problem within the same 
mathematical language of traditional symmetry break- 
ing transitions. Moreover, Wegner introduced a repli- 
cated generating functional F, that can in principle, com- 
pute the exact correlators of the disordered electron gas, 
G(r,r';t) and K(ri, r^, Z\\ r<i, r 2 , z 2 ), provided that the 
non-perturbative aspects of the disorder average in the 
replica limit n — > can be taken correctly Physically, 
the TV-orbital model can be understood as a system of 
metallic grains coupled to its nearest neighbors, forming 
a D-dimensional lattice. The large N limit is equiva- 
lent to the random phase approximation, or a mean-field 
theory of the Anderson model. 

An exact description of the mobility edge has been elu- 
sive, in part, due to the prevailing need to address the 
replica symmetry-breaking substructure of saddlepoint 
manifold which occurs at higher energies away from the 
diffusive regime. Although Wegner's perturbative ap- 
proximation cannot access the transition point, it can 
however confirm particular conductance scaling forms. 
In this paper, an evaluation of F is obtained by exploit- 
ing recent non-perturbative methods developed for spin- 
glass replica field theories with unitary symmetry or bro- 
ken time-reversal invariance @. Previous 0D studies have 



focused on level correlations within the ergodic regime 
uj << 17] (in units of the mean level spacing), however 
the treatment here shall attempt to observe signatures of 
localized states arising from destabilizing replicon spatial 
fluctuations at higher energies. 

Model- Consider the one-particle, single band, tight- 
binding model of Wegner[8j, with N orbitals per site, 
|r, a), a = 1,...,JV on a regular lattice. The pair of 
indices (r, a) will be abbreviated by x, 



(i) 



and defined within a basis that is invariant under the 
local Z 2 gauge transformation, 



|r, a) 



\r,a) 



±1 
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The matrix elements f xx > form a real symmetric 
matrix or RME ensemble, distributed according to 
a normalized Gaussian weight P[fxx>]> with ensem- 
ble averages taken over independent matrix elements: 
M = J D[f]P[f]A[f], assuming {f xx ,\ = and 

\fx\X\' fx2X2 ! \ \^X\X2'^X2X\' ~l~ ^X\X2 ^X\ J X2 f } -^r\— rx' 1 

with disorder width Eq = 4^ r , M r _ r < and K = 1. It is 
necessary to define the normalized inverse, 



4uv 



'/El, such that J2 r > ' 



1. 



Wegner has also introduced the phase invariant ensem- 
ble (PIE), with k = and eigenstates having U(l) gauge 
symmetry Q, 



|r, a) -> e l * r -° \r, a) 



(3) 



The averaged correlators of the RME (PIE) ensemble 
can be expressed as the expectation values of products 
of real (complex) valued spin vectors S(x), with com- 
ponents that are generated by the replica trick S a , 
a, . . . , n, taken in the limit n — > o[j|: 

G(r, r'; z ± ) = [[<r| (z ± - Hy 1 |r')]] = ±i (S a (r)S a (r')) 

(4) 

K(ri ) r' 1) z + ;r 2 ,r' 2 ,z-) = 
(r 1 \(z + -H)- 1 \r 1 ')(r 2 \(z^-H)- 1 \r 2 ')]] (5) 
(S a (ri)S a (ri')S b (r 2 )S b {r 2 ')) 
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with energy and frequency, z± = E(k) ± (uj/2 + id), for 
E(k),u>, 8(z M. In terms of standard transport quanti- 
ties, the mean-free path is given by I = vpT e i, where 
Vpis the Fermi velocity and r e / is the momentum relax- 
ation rate within the Drude picture of electronic scat- 
tering, allowing the disorder strength to be expressed as 
i] = Eq-kN^f) oc l/r e ;, with the density of states at the 
Fermi level A/"(cf). Next, in order to the define a suit- 
able functional integral, the random tight-binding system 
can be reformulated in terms of a continuous, spin-glass 
Hamitonian, 

U sg = TiY] ,S a (r)f r , r ,S a (r')+z±y* S 2 a (r) (6) 

allowing the averages defined through Eq.(|5]) to be ex- 
pressed as, 



(A) = 



fD[S]}A[S}exp(-n sg ) 



(7) 



Wegner Functional- After integrating out / r ,r', a four 
spin interaction is decoupled by a Hubbard-Stratonovich 
transformation and the S a 's are integrated out to obtain 
the celebrated Wegner functional of interacting matrices; 



dQexp (F(Q)) 



F(Q) = --x 



(8) 
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where I is the identity matrix and the trace operation is 
taken in replica space. When the spins of Eq.© interact 
with an infinite-ranged potential, the resulting matrix 
interactions are short-ranged, yielding a OD replica field 
theory which was exactly solved in reference |6| in the 
context of the spherical spin glass problem. However, for 
the iV-orbital model studied here the spherical constraint 
will be supplanted by a saddle point solution for Q at each 
lattice site, or in the continuum limit, an integration path 
in Q-space that must be deformed such that it passes 
through the saddlepoint, which depends on the choice of 
z±. The relevant parts of the solution in ref. Q will be 
reviewed in the next section and the fluctuations about 
the OD saddlepoint will be studied thereafter. 

Exact OD solution- The OD limit, is a single metal- 
lic grain, where the primary complication of concern is 
the non-perturbative evaluation of F in the n — > limit. 
Therefore, the key technical advancement applied here is 
described in the reference jj§], where a Jacobian of Q in 
terms of its eigenvalues Xi, called the Vandermonde de- 
terminant: dQ = ]li<j<fc<ri (^*> _ -\?) G3> auows Eq.© 
to be transformed into a Coulomb gas partition function 
taken in the limit n — > 0: 



Z n = J cxp [H C g(X)] d\i . . . dX n 



(9) 



containing a new Hamiltonian T-L cg describing a one- 
dimensional gas of particles with logarithmic interac- 
tions, often encountered in random matrix theory (lp|. 



H cg {X) = Y J .VTO + /3 5" >|Ai-A 



(10) 



(ii) 



coupling each charge to the single particle potential, 

Fixing /3 — 2 corresponds to the unitary ensemblefioj]; 
the partition function for this system was shown to be 
equivalent to the t iv function of Painleve systems Q, 



n] = — I <lr { 



dec ,f , x 



j=l \<j<k<n 



(12) 



The connection to the OD Wegner functional can be ob- 
tained by taking Xj -> y j y / N/2E$, and t -> zy / N/2E$, 
yielding 



(iv/2i; 2 )("("- 1 )/ 2 +"+^/ 2 / d Vl 
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for which the recurrence relations, for large u = —N/2 
were explicitly solved to yield the following |6j], 
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(14) 

After completing the replica calculation by taking the 
derivative, 

F = — [In Zj = lim (Z n - 1) /n = dZ n /dn\ n=0 . (15) 

the dominant contributions to the functional become, 
F(x) = Nx 

(W) 

where x — z±/Eq. Upon comparison with Eq.([8]), the OD 
saddlepoint Q's in the n — > limit are taken to be 



(x- \Jx 2 - 8) 



= (z ± /E - Jzl/El - 



(17) 
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Accordingly, this agrees with the predicted saddlepoint 
structure for r-independent Q's discussed by Schafer and 
WegnerQ, such that at u) — > 0, Q is bounded by E/E = 

±VE. 

Spatial fluctuations about 0D- The low energy behav- 
ior of the Goldstone or replicon propagator, K(q,w) ~ 
1 / (— iui + Dq 2 ) should be consistent with the well known 
scaling theory of localization^. However, as we will 
demonstrate here, in three-dimensions (3D), at small 
enough E/Eq, the onset of a negative gap should desta- 
bilize the ergodic regime into strong localization, where 
the dc conductivity should vanish. A precise criterion 
where this occurs can be determined by examining dy- 
namical fluctuations about the exact 0D saddlepoint, in 
the n — > limit. Unlike earlier approaches, our func- 
tional is valid at higher energies since the log term of 
F(Q) is treated non-perturbatively, yielding a solution 
very different from the typical Wigner-Dyson correlations 
at lower energies 0, EH- First, consider a fluctuating ma- 
trix field Q'(f) — Q(r) — Qo with coordinates shifted to 
the center of mass frame, r = R + p/2, f = R — p/2 and 
r — f = p. In order to examine the dynamics, the matrix 
field fluctuations can be expanded out to the 2nd order 
in gradients as follows, 



Q'(f) = Q'(R) + fVQ'(i?) + l(£) 2 V 2 Q'(i?) 



Q>(7) = Q'(R) - § VQ'(R) 



l{i) 2 V 2 Q'(R) 



(18) 



The primary contributions to the replicon modes should 
only depend on the non-vanishing 2nd derivative with 
respect to fluctuations, 



SF(Q) = F(Q(r)) - F(Q ) 



p,R 



SF 



SQ'(R)SQ'(R) 



Q'(R) 



(19) 



After taking the appropriate functional derivatives, the 
first term of Eq.® yields, 



sf 1 (q) = -J2Q'(R) 

p,R 



2w p + Avjp 



Q'(R) (20) 



and by defining the mobility constant, l/p p = (p 2 ) 
^2 w p p 2 , the sum over p yields: 
p 



5F 1 (Q) = -'£Q'(R) 



2 + 



Y(R) (21) 



R 



Next, we can perform a similar expansion of the log term 
ofEq.©, 

SF 2 (Q) = E Q'(R) \(E/E - Qo)- 2 ] Q'(R) (22) 



and combining both expressions, the integrated replicon 
propagator becomes: 



K(z±) 



rejr 



d d q 



o p P L q 2 -2 + (z±/E -Q y 



(23) 



It should be noted that q is not the single particle momen- 
tum from the dispersion E(k) — h 2 k 2 /2m — tp, rather 
it is Fourier component of the center of mass coordinate, 
R, of which must be summed over the lattice or in the 
continuum limit, integrated over space. The full con- 
ductivity a should be obtained from summing over both 
advanced and retarded poles, z±, however since K{E) is 
sufficiently long ranged, the short-ranged contributions 
arising from the s summation can be neglected 12]. Tak- 
ing the Fourier transform of integrand in Eq. (|23[) . 



S(R) 



1 ^ f e~ l t R 
4^^ / p- 1 q 2 -2+(z s /E -Qor 2 



s=± 

the conductivity becomes, 



a(u,E) cx uj 2 / R 2 S(R)d d R 



-d a q 

(24) 

(25) 



Thus, in order to ascertain the bounds on a, as in case 
of weak-localization, the divergence of the integral (|2"3")l 
indicates the onset of localization at critical value of the 
scaling parameter t c = E c /Eq, in the dc limit U) — > 0, 
which is independent of the sum over both saddlcpoints, 
z±. 

The mobility edge- In 3D we have, 

K{E) - n h (,,y 2 A\ / tan-^fcF/A^) A 2 > 
4n {V P A J \ tanh-^fc^/A^) A 2 <0 

(26) 

where the replicon energy gap is defined as, 

A 2 = -(2-((£;/ii; ) 2 -8)- 1 ) (27) 

for which A 2 < 0, within (E/E ) E (-y/8,y/8). There- 
fore, K{E) clearly suffers from a strong negative diver- 
gence in the regime kp / ^/JTp'A > 1, which generates the 
condition for localization, 



kp 



p p (2-((E/E ) 2 -8) ') 



> 1 



(28) 



or the mobility edge trajectory in terms of a single scaling 
parameter E/Eq. There are several quantitative limits of 
interest for which the consistency of condition (|27|) can 
be confirmed. At weaker disorder strength (larger E/Eq), 
K{E) > is simply offset by small, negative corrections. 
In the strong disorder limit Eq — > oo, K(E) acquires a 
large negative gap in a manner similar to the de Almeida- 
Thouless instability of mean-field spin glass systems [Uj]. 
An important distinction here is that I considered spa- 
tial fluctuations about the 0D saddlepoint rather than 
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for 
> 



replica symmetry-breaking fluctuations about the replica 
symmetric saddlepoint. 

In order to make a connection with known results of 
the single orbital Anderson model with unitary symme- 
try, the scaled critical width was determined numerically 
as Efi sa 18.3[l4T], In the TV-orbital model, the mobility 
inverse l//i p can be re-expressed in terms of transport 
parameters by dimensional analysis: l//i p oc AEq 
some positive constant A, and Eq, follows from E 
in Eq.dHJ), yielding {E^} 2 = 15/(8,4). 

Discussion- My primary focus here was to demon- 
strate critical localization within the framework of a OD 
replica field theory beyond the low energy limit. Conse- 
quently, the precise determination of the transition point 
was not possible with earlier low energy, non-linear sigma 
model approaches that are only valid within the diffusive 
regime [1 If. Hence, I have presented the gradient expan- 
sion of spatial fluctuations about the single grain saddle- 
point as a new mean-field, hydrodynamic approach which 
can access the breakdown of the ergodic regime into local- 
ized states within a mobility edge trajectory. These spa- 
tial correlations can be quantitatively bounded within a 
Ginzburg criterion, where the variance of the spatial fluc- 
tuations should be compared to the square of its average 
value. The negative divergence of the integrated replicon 
propagator clearly demonstrates a breakdown of the dc 
conductivity and the fundamental necessity of restoring 
non-mean-field corrections, which are beyond the scope 
of the gradient expansion. 

The single grain saddlepoint utilized here is a math- 
ematically exact starting point for studying Anderson 
localization, facilitating investigations of complexities 
severely lacking in conceptual understanding such as the 
replica symmetry-breaking substructure of the localized 



regime 15], which is common to both spin-glass and local- 



ization phenomena. Hence, the methods introduced here 
can be further refined to include fluctuations in both real 
and replica spaces by taking the expression of the OD sad- 
dlepoint solution T IV [n] given by Eq. (fT4|) . which precedes 
the n — > limit, and including replica symmetry breaking 
dynamics, thereby developing a stability criterion along 
the lines of de Almeida and Thouless[13j. Thus, for a 
given symmetry class of disordered systems, a different 
OD dQ Jacobian will lead to a different type of Coulomb 
gas ensemble, which can in principle can be evaluated in 
terms of random matrix averages. 

Recently, Anderson localization has played a central 
role in determining the robustness of surface states of 
topological insulators [16j. It particular, it has been dif- 
ficult to reconcile the topological terms that represent 
magnetic field or Berry curvature effects [13] within the 
sigma model approach and access higher energy localized 
states. Beyond the approximation of a spatially uniform 
saddlepoint, the model studied here at f3 — 2 is of the 
unitary symmetry class, which is the simplest case to con- 



sider because it avoids the complexities of Pfaffians that 
arise in the orthogonal ensemble f3 — 1. Thus, an exten- 
sion of the methods here may offer an alternative route 
towards a solution of the the elusive the integer quantum 
Hall plateau transition, which is a mobility edge problem 
in the presence of a magnetic field through a modifica- 
tion of F to include a topological term. In that case the 
2D quantum hall system is delocalized at a critical value 
of the topological winding number, = 7r in the diffu- 
sive regime, however a more complete description would 
involve localization criterion that includes both and 
the scaling parameter E/Eq into an inequality similar to 
Eq.([23). 
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